Ferromagnetism in the Fe-substituted spinel semiconductor ZnGa2 04 
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Motivated by the recent experimental observation of long range ferromagnetic order at a relatively 
high temperature of 200K in the Fe-doped ZnGa2 04 semiconducting spinel, we propose a possible 
mechanism for the observed ferromagnetism in this system. We show, supported by band structure 
calculations, how a model similar to the double exchange model can be written down for this system 
and calculate the ground state phase diagram for the two cases where Fe is doped either at the 
tetrahedral position or at the octahedral position. We find that in both cases such a model can 
account for a stable ferromagnetic phase in a wide range of parameter space. We also argue that 
in the limit of high Fe^^ concentration at the tetrahedral positions a description in terms of a two 
band model is essential. The two eg orbitals and the hopping between them play a crucial role in 
stabilizing the ferromagnetic phase in this limit. The case when Fe is doped simultaneously at both 
the tetrahedral and the octahedral position is also discussed. 
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I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) are recently 
being intensively studied in connection with their possi- 
ble application for spintronic devices Jj. Special atten- 
tion has been devoted to the III-V semiconductors 0, Q 
which develop long range ferromagnetic order with Curie 
temperatures of about lOOK upon doping with a low con- 
centration of magnetic impurities like Mn. Since spin- 
tronic applications would become widely accessible if the 
ferromagnetism is achieved at room temperature, there 
is a continuous search for new materials with high Curie 
temperatures. 

In a very recent experiment Risbud and coauthors 
tried to dope Fe into ZnGa204 by preparing a solid so- 
lution [ZnGa204 ]i_a;[Fe304]a; of ZnGa204 and Fe304 
with x=0.05,0.10 and 0.15. Long-range magnetic order 
was observed in all the three samples with Curie tem- 
peratures up to 200K as well as ferromagnetic hysteresis 
of the magnetization at low temperatures. Interestingly, 
the saturation magnetic moment which should be ifiB 
per Fe304 unit, is about IfiB instead. This has been 
interpreted as an indication of certain fraction of Fe 
not contributing to the ferromagnetic long range order 
(LRO). These authors have also performed Mossbauer 
experiments in order to ascertain the oxidation states of 
Fe in the host semiconductor ZnGa204. For a, x = 0.15 
doped sample they observed the presence of only Fe'^+ 
states with some of them displaying a paramagnetic sig- 
nal (doublet) and the rest showing magnetically ordered 
Fe'^"'". We note at this point that in the context of Fe304, 
it has already been argued that above the Verwey tran- 
sition ^ an Fe^+ cation can be viewed as an Fe^+ ion 
plus a delocalized electron. We deliberate further on this 
point in section VII. 

Motivated by these observations and the difference 
between this system and the III-V semiconductors, 
we investigate in what follows the underlying mecha- 



nism of the long range ordering of Fe ions in doped 
ZnGa2 04 considering various possible limits of the prob- 
lem in terms of an effective model which is based on 
the band structure calculation of the system. The band 
structure calculation, which will be discussed in the next 
section, gives the information about the position and na- 
ture of the Fe bands, their hybridization with the bands 
of the parent compound as well as the active bands at 
the Fermi level which will help us to construct an ef- 
fective Hamiltonian for the system. The ground state 
magnetic phase diagrams of this effective Hamiltonian 
are then calculated using the parameters derived from 
the band structure calculation. 

The host semiconductor ZnGa2 04 has a spinel crys- 
tal structure AB2O4 with two cation sites: Zn^+ (A) in 
a tetrahedral co-ordination and Ga'^+ (B) in an octahe- 
dral co-ordination of oxygens. Fe3 04 has, on the other 
hand, an inverse spinel structure with a chemical compo- 
sition Fe^+[Fe2+,Fe^+]B04^. When Fe is substituted in 
ZnGa204 via the solid solution [ZnGa204 ]i^x[Fe304]x, 
it can either replace Zn in the tetrahedral position or Ga 
in the octahedral position or both. Here we will con- 
sider the following two cases: i) all substituted Fe are 
in tetrahedral positions and ii) all substituted Fe are in 
octahedral positions. We also assume the most general 
case, namely, that Fe can have both Fe'^"'" and Fe^"*" ox- 
idation states irrespective of whether it is in tetrahedral 
or octahedral position and that it is always in a high- 
spin state with spin=5/2 and spin=2 respectively We 
will briefly outline the case where Fe ions are in both 
tetrahedral and octahedral positions at the end. 

The paper is organized as follows. In Sec. II we discuss 
our band structure calculations and present the density of 
states for the case of Fe doped into tetrahedral positions 
in ZnGa204 . Out of these results we get insight about 
the active orbitals at the Fermi level which are relevant 
for our model Hamiltonian. In Sec. Ill we investigate 
three different limits of the case with Fe occupying tetra- 
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hedral positions and motivate a model similar to double 
exchange for this system. In Sec. IV we present the 
model Hamiltonian and calculate the magnetic phase di- 
agram. Sec. V deals with the effect of the Coulomb 
correlation on the phase diagram results. In Sec. VI we 
study the case of Fe doping in the octahedral positions 
and finally in the last section we discuss our results and 
make a comparison with the experimental observations. 



II. BANDSTRUCTURE CALCULATIONS 

ZnGa204 is a direct band gap semiconductor with an 
energy gap of about 4.1eV |a, Q- Previous band struc- 
ture calculations for ZnGa2 04 showed that the va- 
lence states right below the Fermi level are mostly of 
oxygen character with the contribution of Zn and Ga be- 
ing very small. In order to investigate the effect of the 
doping of Fe in the band structure of ZnGa204 , we con- 
sidered modified unit cells of ZnGa204 with different Fe 
content, i.e. Fe substituting Zn in tetrahedral sites in an 
Fe:Zn ratio 1:2 (50% doping of Fe) and Fe substituting 
Ga in octahedral sites in an Fe:Ga ratio 1:4 (25% doping 
of Fe). These percentage doping of Fe in Zn (or A) sites 
and in Ga (or B) sites are close to the experimental ones 
with X — 0.15 in [ZnGa204 ]i_a;[Fe304]2; formula unit 
which are 45% when Fe substitutes Zn only and 22.5 % 
when Fe substitutes Ga only. Here we will illustrate the 
tetrahedral substitution. 

ZnGa204 crystallizes in a normal spinel structure with 
space group Fd3m and the primitive (rhombohcdral) 
unit cell contains two formula units. In order to sub- 
stitute Fe in one of the Zn positions, one has to make 
the two Zn positions in the unit cell nonequivalent. One 
maximal subgroup of Fd3m which allows for this sub- 
stitution is FASm. In this new space group we have in 
addition to Zn, Fe and Ga, two nonequivalent oxygen 
positions Oi and O2. Fe and Zn are surrounded by Oi 
and O2 ions respectively in a tetrahedral environment, 
whereas Ga is in an octahedral surrounding with three 
Oi and three O2 atoms '§]. 

We have performed ab initio density functional theory 
calculations for Fe doped ZnGa204 in the F43m symme- 
try within the local spin density approximation (LSD A) 
using the linearized augmented plane waves (LAPW) as 
basis set 0. In Fig. we present the spin polarized 
density of states (DOS) for Fe doped ZnGa204 where Fe 
substitutes 50% of tetrahedral Zn atoms in the unit cell 
as described above. Total density of states of majority 
(up) and minority (down) spin states for all the atoms 
are given in Fig. ^a) while in Fig. ^b) we show the 
partial density of states for Fe rf— states projected into 
t2g and Bg symmetries in both spin directions. We ob- 
serve that the Fe d-states appear to be mainly located in 
the band gap of the host semiconductor together with a 
non-negligible contribution of Oi p-states while Zn, Ga 
and O2 contributions are well down into the valence band 
and are negligible at the Fermi level (Fig. ^a)). The Fe 
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FIG. 1: (Color) Majority (up) and minority (down) spin 
density of states for Fe doped into the tetrahedral Zn position 
in ZnGa2 04 . (a) shows the total density of states and (b) the 
partial Fe-d density of states as explained in the text. Here 
01 and 02 denote two nonequivalent oxygen atoms of the 
modified unit cell. 



DOS shows the expected eg-t2g splitting for a transition 
metal ion in a tetrahedral crystal field. The majority 
spin (spin up) Cg and t2g states are completely filled and 
appear far below the Fermi level whereas the minority 
spin (spin down) Cg states are partially filled and are at 
the Fermi level (Fig. ^b)). It is important to note that 
the spin exchange splitting 2.6eV) is much larger than 
the crystal field splitting (~ 0.5eV) in this system. These 
observations of the band structure calculation are crucial 
to build up an effective model for this system, as we will 
see below. 



III. MODEL : Fe IN TETRAHEDRAL POSITION 

Based on the electronic structure calculations des- 
cribed above for the doped system, we now motivate a 
possible mechanism for the ferromagnetic LRO observed 
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(a) (b) 

FIG. 2: Schematic energy level diagram showing the 3d states 
of doped Fe in the band gap of ZnGa2 04 semiconductor, (a) 
Fe^^ and (b) Fe^^ in a tetrahedral crystal field. V.B. and 
C.B. denote valence band and conduction band of the semi- 
conductor respectively. 



in Fe doped ZnGa204 . We note at this point that the 
Fe concentration in this system is not very low as for 
x=0.15 in the [ZnGa204 ]i_a;[Fe304]2: formula unit the 
Fe:Zn ratio calculated from the nominal composition is 
1:1.9 (3x:l-x) or close to 1:2 though Risbud et. al. re- 
ported to have observed a Fe:Zn ratio of 1:3 Hence 
the system is not likely to be a candidate for description 
in terms of a Kondo impurity model |llj| . It is rather 
much more reasonable to assume that the electrons hop 
between doped sites via mainly the oxygen orbitals in the 
same way as in the double exchange mechanism in the 
manganites |il2]. We would like to note here that in the 
case of manganites each Mn ion is surrounded by oxygen 
atoms in an octahedral environment and these octahe- 
dra are corner shared. The double exchange mechanism 
then operates through a path of the type Mn-O-Mn giv- 
ing rise to a ferromagnetic order among the Mn spins. 
In the present system we observe that for the doping 
considered in our band structure calculation there ex- 
ist paths of the type Fe-Oi-Oi-Fe where Oi is one of the 
two nonequivalent oxygens which surrounds Fe in a tetra- 
hedral environment as described in the previous section. 
This assumption of electron transport via oxygen orbitals 
is supported by our band structure calculations where we 
see that there exists a small but finite hybridization be- 
tween the doped Fe and Oi at the Fermi level (see Fig. 
^a)). Note that the contribution of Zn, Ga and O2 is 
negligible at the Fermi level and hence are very unlikely 
to take part in the electronic transport. 

A double exchange like mechanism mediating fer- 
romagnetism has been proposed and is being seri- 
ously investigated lately for various diluted magnetic 
semiconductors [13I IT^ ITs f . First principles calculations 
have also been seen to support such a mechanism for 



ferromagnetic order in some of these systems like ZnO 
based DMS, Ga(Mn)As [Hill. Indeed, without such a 
long range transport of electrons, the ferromagnetic LRO 
(Long Range Order) observed in these systems would be 
difficult to account for, a view shared in other theoretical 
analyses^ of LRO in DMS. We would like to emphasize 
though that the Fe doped ZnGa204 that we are consid- 
ering here is not strictly a DMS as the doping level is 
quite high. 

Let us consider first the case where all the doped Fe are 
in tetrahedral positions. In Fig. 2 we draw a schematic 
energy level diagram of Fe id orbitals together with the 
valence and conduction band of the host semiconductor. 
We discuss in the following three limits of the problem, 
(i) all or most of the Fe are Fe^+, (ii) Fe is in both Fe^+ 
and Fe'^+ oxidation states and (iii) the limit of high Fc^"*" 
concentration. 



A. High Fe concentration limit 

In the situation when neighbouring doped tetrahedral 
sites are all Fe'^"'" , there are 5 electrons of the same spin 
in each site due to Hund's rule. If the spins at these 
neighbouring doped sites are ferromagnetically aligned 
then hopping of an electron between sites is blocked by 
the Pauli principle. If instead, they are antiferromagnet- 
ically aligned the system gains superexchange energy by 
a virtual process of electron transfer between the Fe ions. 
Therefore, when all Fe ions are in a Fe'^"'" state, an anti- 
ferromagnetic (AFM) alignment of spins is energetically 
preferred. 



B. Mixed Fe^+ and Fe^+ 

Let us examine now the case where both Fe^^ and 
Fe'^"'' are present in neighbouring sites (see Fig.O. If the 
5 electrons in Fe^"*" (all aligned) are in an antiferromag- 
netic configuration with the spins of the neighbouring 
Fc'^+ (Fig-El^)), then an electron in Fc^+ has to pay an 
amount of energy equal to the Hund's coupling {Jh) in 
order to hop from one Fe^+ site to a neighbouring Fc'^^ 
site. In the limit of large Jh this is practically forbid- 
den and the system will try to gain the superexchange 
energy, approximately -j^, where t is the appropriate 
hopping integral between the relevant orbitals On 
the contrary, in a ferromagnetic arrangement(Fig.|3|^b)), 
the minority spin electron in Fe^^ can move to a neigh- 
bouring Fe'^+ site without paying extra energy and the 
system gains kinetic energy (KE) in this process. From 
the above discussion, it is evident that when Jh is large, 
the system will prefer to be in the ferromagnetic state 
rather than the antiferromagnetic one. If Jh is moder- 
ate then all these energy scales are comparable and the 
competition between kinetic energy, superexchange en- 
ergy (SE) and Jh will decide the phase boundaries. 
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Fe2+ Fe 
(b) 

FIG. 3: Electron hopping between Fe^""" and Fe"'"'" at neigh- 
bouring tetrahedral doped sites with (a) antiferromagnetic 
ahgnment of spins (b) ferromagnetic alignment of spins. 



It is useful to make a note at this point that the Fe^"*" 
ion with 3 electrons in its eg orbitals is likely to be Jahn- 
Teller (JT) active, i.e. the doubly-occupied orbital stabi- 
lizes over the singly occupied one, whereas Fe^+ with two 
electrons in the two Cg orbitals is not JT active. In this 
case one can then work with only one eg orbital for Fe^"*" 
(the Jahn- Teller stabilized one) with electrons hopping 
through this orbital. The mechanism of magnetic ex- 
change due to electron delocalization goes through with- 
out loss of generality as outlined above (with possible re- 
duction in the overall eg bandwidth, which can be scaled 
away) . The third eg electron in the Fe"'"^ site and the cor- 
responding Cg orbital can be ignored. But as we will see 
in the following, this single orbital model is not sufficient 
to describe the case where all (or most of the) doped Fe 
are in a Fe^+ state. One has then to take into account 
both eg orbitals and the hopping among them. 



C. High Fe concentration limit 



Here we consider two neighbouring tetrahedral Fe ions, 
both in Fe^"*" configuration. Since each of them have 3 
electrons in their eg orbitals, one eg orbital is full and the 
other has one electron. Therefore the only way electrons 
can hop is via this half-filled orbital. In this case, there 
are several possibilities arising from the relative values of 
the JT stabilization energy (A) and the bandwidth {W) 
of the Cg bands (A << Jh for the system under consid- 
eration). Consider the situation when A > W^. In this 
case ferromagnetism is inhibited if the hopping matrix is 
diagonal. If, however, there exists off-diagonal hopping 
t^^ (which, in general, depends on the orbitals involved 
and the symmetry of the lattice), then ferromagnetism 
could stabilize via a virtual hopping with a gain of fer- 
romagnetic exchange energy (t^^)^/A. This FM phase is 
not driven by KE as in the double exchange mechanism. 
However, there is also a competing AFM phase that gains 
superexchange energy of order / (A + Jh)- 

In the limit A < W ^ the FM state is driven by the 
double exchange mechanism as in the Fe^+ - Fe"^"*" mixed 
configuration. The FM state is stabilized by the KE of 
the Qg electrons since the superexchange energies are less 
than the KE. The underlying ground state, though, will 
be different when only diagonal hopping is allowed. In 
this case, in order to gain the KE, the Fe^+ ions will 
remain in a cooperative, staggered JT distorted arrange- 
ment which costs additional energy which depends on A 
and may not be stable if AFM superexchange energy is 
larger [l^ . 

In Fig. ^ we show the situation arising in the case 
A = with the eg orbitals in a cubic environment and 
nonzero overlap among them. One electron (either up 
or down) from the doubly occupied eg orbital of Fe^+ on 
one site can always move to the singly occupied eg orbital 
on the next Fe^"*" site. In case of antiferromagnetic align- 
ment (shown in Fig. EJa)), the hopping of an electron 
from one Fe^"*" site to the next costs Jh amount of energy. 
The superexchange mechanism through virtual hopping 
is the only energy gain. In the ferromagnetic state (Fig. 
01b)), however, the initial and final states are degenerate 
and the system gains kinetic energy due to resonance. It 
is also evident that the physics for large Fe^"*" concentra- 
tion is different from the low Fe^+ concentration limit. 
In the low concentration limit, an effective single orbital 
model captures the physical situation well, while in the 
high concentration limit, two bands are crucial for its 
understanding. This situation is somewhat reminiscent 
of the manganites where in the electron-doped side (hole 
concentration x > 0.5) the two eg orbitals and the hop- 
ping between them play a crucial role in determining |20j 
the competition between the different magnetic phases 
whereas in the hole-dope side (x < 0.5), a model with 
only the Jahn- Teller stabilized single orbital is adequate. 
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FIG. 4: Electron hopping in (a) antiferromagnetic phase (b) 
ferromagnetic phase when both neighbouring doped tetrahe- 
dral sites are Fe^^. 



IV. HAMILTONIAN 

Based on the above discussion we write down 
the following effective Hamiltonian for Fe-doped 
ZnGa204 assuming a cubic environment 

H = Hq + Hint 
^0= ttf^aaCjPa-JH^S^.a, 



+ JaF ^ Si.Sj 



<ij> 



Hint = U^^n^airiiai +U' ^ niariip 



(1) 



Here we treat the t2g electrons as localized and eg as 
itinerant because the density of states (Fig. 1) clearly 
shows that the Cg down band is at the Fermi level, 
whereas the fully filled t2g up band is well below. Elec- 
tronic transport, therefore, involves the electrons in the 



Cg band primarily. The t2g bands are well removed 
from the Fermi level and the eg-t2g overlap is negligi- 
bly smalljl^. The t2g electrons, under these conditions, 
provide a localized magnetic background, to which the 
itinerant eg electrons are coupled through the Hund's 
exchange. 

The first term in Hq describes the kinetic energy with 
t"j^ being the anisotropic hopping integral between two eg 
orbitals[22|2j|3|- Here i,j are site indices and a, (3 =1,2 
are eg orbital indices. The second term is the Hund's cou- 
pling term between the localized i2g spins and the itin- 
erant eg spins and the last term represents the antiferro- 
magnetic superexchange coupling between neighbouring 
t2g spins. The first and second term in Hint define the 
onsite intra- and inter-orbital Coulomb repulsion with 
U and U' being the corresponding interaction strengths. 
In the half-filled situation when we have all Fe in Fe'^"'" 
states the ferromagnetic phase is blocked by the Pauli 
principle and an antiferromagnetic phase is favoured as 
we discussed above. The third term in Hq representing 
the antiferromagnetic superexchange satisfies this limit. 

First we discuss the model without considering the 
Coulomb correlation terms given by Hint (Eq. The 
correlation and their effects will be discussed later in 
detail. We treat the t2g spin subsystem with a mag- 
nitude 3/2 semiclassically as it is also the standard prac- 
tice in the case of manganites We do not, how- 
ever, make the further assumption prevalent in mangan- 
ite and DMS literature i.e., — > oo. Indeed, such 
an assumption would preclude the presence of a Fe"*"^ 
state. Jh in the foregoing is treated as a parameter 
and its value, as gleaned from the spin splitting ob- 
served in the band structure calculations, is typically 
large (about 2.(!>eV). Assuming an uncanted homoge- 
neous ground state, we choose Si = Soexp{iQ.ri) where 
So — 3/2 and Q = (0,0,0) for the ferromagnetic phase 
and Q = (tt, tt, tt) for the antiferromagnetic phase. With 
this choice the first two terms of the Hamiltonian Hq (1) 
reduce to 



a/3 t 



22 



JHSo^cl^^Ck+Qat 
+ JH5o^4^^Cfc+QaJ. (2) 

-2t{coskx + cosky) 



f21 



V3 



t{coskx — cosky) 



-t(coskx -\- cosky) 



-tcosk^ 



(3) 



Here 1 corresponds to dx2-y2 and 2 to (12^2 orbital 
and t is the magnitude of the hopping integral between 
two neighbouring d.j;2_y2 orbitals in the x,y direction. 
The superexchange contribution to the Hamiltonian is 
given by 

EsE = ^^^{2cose,y + cose,) (4) 
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FIG. 5; Ground state phase diagram of the Hamiltonian {U = 
U' = 0) in the y — JHSo/t plane where y is the electron 
concentration and with JafSq/I = 0.05. 



FIG. 6: Phase diagram for Fe in a tetrahedral environment in 
the y — JafSq/I plane considering both eg orbitals and the 
hopping among them in the model in the limit U — U' = 0. 



where 9xy and 9^ are the angles between neighbouring 
spins in the xy plane and in the z-direction respectively. 
9xy — Oz = for the ferromagnetic phase and 6xy — Qz ^ 
TT in the antiferromagnetic phase. These two angles could 
be different from tt or in general and allow for canting. 

We diagonalize the Hamiltonian (Eq.|2J) at each k point 
on a finite momentum grid and calculate the ground state 
energy for ferromagnetic and antiferromagnetic states in 
their uncanted spin configurations. The magnetic struc- 
ture with minimum ground state energy is determined 
for each set of parameters (y, Jh and Jaf), the two e^ 
orbitals are taken to be degenerate presently (A = 0). 
Here y is the tg electron concentration, y = 0.5 corre- 
sponds to the limit where all Fe ions are in their Fe'^+ 
state and y — 0.75 corresponds to all Fe in Fe^+ state. 
In Fig. [3 we show the ground state phase diagram in the 
y — JHSo/t plane with JafSq/I — 0.05 as an illustra- 
tion. However, this value of JApS^/t is varied in a wide 
range to obtain the phase diagram given in Fig. |3 In the 
above figure we see that at y=0.5 where all Fe ions are 
in their Fe''+ state the system is antiferromagnetic at all 
values of the Hund's coupling Jh as we expected since 
the ferromagnetic state is blocked by the Pauli exclusion 
principle in this limit. As we increase the concentration 
of Fe^+ , the electron concentration increases in the down 
spin band which can hop from site to site and the system 
gains kinetic energy. Due to the competition, modulated 
by the value of J/f, between the kinetic energy which 
favours a FM configuration and superexchange energy 
which favours an AFM state, a ferromagnetic phase is 
indeed stabilized over the antiferromagnetic one for mod- 
erate to high concentration of Fe^+ . As the value of Jh 
is increased, the ferromagnetic phase becomes broader 
and at very large Jh the ferromagnetic region becomes 
almost independent of Jh ■ 



In Fig. El we present the ground state phase diagram in 
the y — JAFSo/t plane at a typical value of JjjS'o/t = 25.0 
which is again estimated from the electronic structure 
calculation described above. Note that there exists a 
wide region in parameter space where the ferromagnetic 
phase is stabilized. 

In the limit of non-degenerate e^ orbitals, we examine 
the situation for > A > W. As argued earlier, there 
exists the possibility of a ferromagnetic phase via dou- 
ble exchange here too in the region of mixed Fe+^-Fe^"^ 
shown in Fig. 7. In this case, the AFM state reappears as 
the Fe"*"^ concentration increases because of reduced ef- 
fective hopping. The phase diagram is symmetric about 
y = 0.625 and the regions of stability of these phases 
are nearly independent of A for A > W a.s expected. 
The rather interesting possibilities involving orbital or- 
der have not been discussed here. The orbital order can 
be generated by the anisotropic hopping as well as the 
JT distortion. It can also be enhanced by the Coulomb 
correlations poj. 



V. COULOMB INTERACTION 

The onsite intra- and inter-orbital Coulomb interaction 
terms given by Hint (in Eq. ^ are treated in the mean- 
field theory. Neglecting the fluctuation effects, we write 

nilahi2a' =< flla > rti2<T' + < n2a' > flu^- < flia X 

n2(7' >, the last term preventing double counting. The 
averages < ni| >, < nij >, < n2-f >, < fii^ > are calcu- 
lated from the eigenvectors iteratively through successive 
diagonalization of the Hamiltonian. Self-consistency has 
been achieved when all the averages < hi^a-,a > and the 
ground state energy converge to within 0.01 % or less. 

It is well known [13.12^ that in the large Jh limit, the 
Coulomb repulsion U between up and down spin elec- 
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FIG. 7: The ground state phase diagram for Fe in a tetrahe- 
dral environment in the y — A/t plane at a typical value of 
JuSo/t = 25.0 and JafSI/I = 0.05 



FIG. 8: Phase diagram of the Hamiltonian for Fe in a tetra- 
hedral environment showing the effect of U' (see Eq. on 
the ferromagnetic and antiferromagnetic phases 



VI. Fe IN OCTAHEDRAL POSITION 



irons at the same orbital is ineffective in the mean-field 
theory. In this limit doubly occupied orbitals are energet- 
ically costly and generally avoided. However, in the case 
of Fe"*"^ , one of the Cg orbitals has to be doubly occupied. 
In that case, in the mean- field type argument, a replace- 
ment of by J/f -f [/ in the doubly occupied orbitals 
takes care|26l| of this repulsion. One then absorbs U in 
the value of Jjy appropriately and does not consider it 
explicitly. In the following, therefore, we keep the value 
of U to be zero. There is, however, a strong effect of 
the inter-orbital Coulomb interaction ([/') on the phase 
diagram as we discuss in what follows. 

In Fig. 8 we present the effect of inter-orbital Coulomb 
interaction {U') on the ground state magnetic phase di- 
agram for typical values of JhSq/I and JafS'^/I, keep- 
ing A = 0. As y increases, the FM phase appears as 
in the previous figures. At large U' and when almost 
all the Fe ions are in the +2 valence state, the AFM 
phase reappears at the right top corner of Fig. 8. In 
the presence of inter-orbital Coulomb interaction the en- 
ergies in the high Fe^^ region are primarily dominated 
by localized exchange interactions. The competing in- 
teractions now have the energy scales —jp- and — 
In the limit U' > Jh (Jh + U, ii U is considered ex- 
plicitly), the second term would provide extra gain in 
energy and the AFM phase should stabilize. A transi- 
tion from FM-^AFM will therefore occur as U' exceeds 
J// for y = 0.75. As y reduces from 0.75, a larger U' is 
required for the transition leading to the region of AFM 
at the top right hand corner in Fig. 8 as shown. 



So far we have investigated the case where the doped Fe 
ions are in the tetrahedral Zn positions of ZnGa2 04 with 
Fe^"'' and Fe^+ valence states. Let us now examine the 
case when Fe is doped into octahedral Ga positions. In 
the octahedral crystal field the energy levels of Fe will 
be split into a triply degenerate set of t2g levels and a 
doubly degenerate set of eg levels. The t2g levels in this 
case have lower energy than the eg levels, contrary to 
the tetrahedral case. With this arrangement of orbitals, 
the extra (6th.) electron in Fe^+ will occupy the t2g 
level. The preliminirary band structure results with 25% 
Fe substitution in the Ga sites corroborate this sce- 
nario. The overlap integrals between the t2g orbitals are 
calculated as usual from the Slater-Koster integrals [2^. 
In this scenario, we consider the t2g electron as itinerant 
for reasons similar to the ones discussed in the case of 
Fe in tetrahedral position. We, therefore, use the overlap 
integrals between the t2g orbitals in the kinetic energy 
term of Hamiltonian Eq. |21 We note at this point that 
in the spinel crystal structure of ZnGa204 these octa- 
hedral centers, occupied by the Ga atoms, are arranged 
in a tetrahedral fashion among themselves and hence are 
geometrically frustrated. We have not considered this ge- 
ometrical frustration in the present calculation because 
the experimentally observed ratio of Fe to Ga is fairly 
small (1:6). The Fe atoms are assumed to be arranged 
in a cubic environment for the present calculation and 
hence there is no frustration. 

In the kinetic energy term of the Hamiltonian given by 
Eq. Qthe orbital indices a and /3 now take the values 1,2 
and 3 which represent xy,yz and zx orbitals respectively. 
This Hamiltonian will reduce in k-space to the form of 
Eq. 121 as in the case of tetrahedral doping except from 
the fact that now e'^^ is a 3x3 matrix with the elements 
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FIG. 9: Phase diagram for the ground state of the Hamilto- 
nian (with U = U' — 0) for Fe in an octahedral environment 
in the y — JhSq/I' plane with JafSq/I' = 0.05. In the inset 
we show the shift of the AFM-FM phase boundary compared 
to that of the tetrahedral case (x-axis is rescaled such that 
y' = corresponds to all Fe^"*" and y' — Ito all Fe^"^ in order 
to make a comparison) 



given by 



,22 



,33 



—2t' [coskx + cosky) 
—2t' [cosky + coskz) 
—2t' [coskx + coskz) 



(5) 



Here t' is the magnitude of the hopping integral be- 
tween the neighbouring 7r-bonded xy orbitals in the x,y 
direction. Note that the inter-orbital overlaps turn out 
to be zero in this case Following similar procedures 
outlined in the tetrahedral case we calculate the ground 
state magnetic phase diagram and observe that (as shown 
below) the ferromagnetic phase is stable in a wide range 
of parameter space in this case as well. However, we have 
not considered the Coulomb interactions {Hint in Eq. ^ 
in the present calculation of the phase diagram. 

In Fig. 1^1 we present the ground state phase diagram in 
the plane y — JuSo/t' with a fixed value of JafSq/I' — 
0.05. The three t2g orbitals have been taken as degen- 
erate (i.e., A = 0) and y is the electron concentration 
in the t2g levels ranging from y = 1/2 (corresponding to 
all Fe3+) to y = 2/3 (all ¥e^+ , 4 electrons in t2g). A 
comparison with the earlier phase diagram, for Fe in a 
tetrahedral position (Fig. [S]), is shown in the inset. We 
observe a shift of the AFM-FM phase boundary towards 
higher concentration of Fe^+ in the octahedral case. The 
absence of off-diagonal hopping among the t2g orbitals 
(Eq. |SJ) reduces the effective KE gain in the double ex- 
change mechanism. This, in turn, makes the FM phase 
less stable compared to the tetrahedral case. The an- 
tiferromagnetic phase, therefore, stabilizes over a wider 
region in the phase diagram. 



Nevertheless, we still observe a stable ferromagnetic 
phase in a moderate to high range of doping by Fe^"*" 
ions. Finite JT splitting of i2g orbitals may have inter- 
esting effects on the stability of the ferromagnetic phase 
in the high Fe^+ limit. Since the inter-orbital hopping is 
zero, even a small but finite JT splitting may induce a co- 
operative staggered JT distorted order of the Fe^"*" ions in 
the ferromagnetic phase to maximize the KE. This would 
make the ferromagnetic phase increasingly destabilized 
against the AFM phase, depending on the value of A, as 
well as against defects and other impurities present in the 
system jl9|. Coulomb correlations, particularly the inter- 
orbital Coulomb correlation ([/'), is expected to have a 
strong effect in stabilizing the AFM phase in this case 
due to the reduced mobility of the carriers in the t2g or- 
bitals as argued above. Following the arguments in Sec. 
V, the AFM phase at the right top corner in Fig. |S1 is 
expected to appear at a lower value of U' now. 



VII. DISCUSSION 

We observe from the above study of Fe doped 
ZnGa204 that ferromagnetism could result from a deli- 
cate competition between double exchange, favouring the 
FM spin order, and superexchange, favouring AFM or- 
der, when both Fe^"*" and Fe^+ valence states are present. 
Though the present experimental situation in this sys- 
tem does not make the limit of all Fe^^ in tetrahedral 
positions a very relevant one, there is nevertheless a 
rich underlying physics connected with the interplay be- 
tween dominant Jahn- Teller, double exchange and su- 
perexchange interactions in this limit. Finally, a special 
note is in order for the case where the doped Fe ions go to 
both tetrahedral and octahedral environments replacing 
Zn and Ga ions respectively. In this mixed situation, it 
is necessary to first find out the extended and localized 
states from a careful density functional calculation. As 
the eg-t2g orbital overlap is negligible, it could be pos- 
sible that in the mixed situation the double exchange 
mechanism would operate within the Fe ions that belong 
entirely to one kind of crystalline environment (i.e., ei- 
ther in the tetrahedral or in the octahedral positions) and 
hence only a fraction of the Fe that are doped into the 
system would take part in the ferromagnetic long range 
order However, the model is expected to work as 
long as there are mobile electrons coupled to a relatively 
localized spin background. 

We have been discussing in this work the carrier- 
mediated ferromagnetism which is believed to explain the 
magnetic properties of various dilute magnetic semicon- 
ductors with available free charge carriers 

MM For 

relatively localized systems with no free carriers other 
alternative mechanisms such as that of the bound rnag- 
netic polaron model have also been proposed |23, l28j . 
One such example is Ga(Mn)N where it has been sug- 
gested that the Mn ions are in d^ configuration plus a 
localized hole and this localized hole forms a singlet with 
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a Mn d-electron (Zhang-Rice polaron) which then moves 
through the Mn sublattice and mediates ferromagnetic 
order j29ll3(il|. 

The system we consider in this work is different from 
other III-IV or II- VI semiconductors in some respects. 
First of all the doping concentration (50% or 25%) is 
quite high compared to, for example, the 6% doping in 
Ga(Mn)N. Secondly there is a possibility of Fe having a 
mixed valence state of Fe^"'"/Fe'^"'". Thirdly the existence 
of two kinds of crystalline environment to which Fe can be 
doped is not present in the known DMS semiconductors. 

DMS Systems like Ga(Mn)As have indeed shown de- 
pendence of Tc on the so called antisite defects. Long 
range order in DMS systems is known to show sensitivity 
to disorder as well. But in the present context, disorder 
(which we did not consider) may not be that crucial for 
the underlying mechanism of magnetism proposed. Un- 
like the usual DMS materials, the doping is fairly high to 
be in the impurity dominated regime. The disorder may 
affect the values of exchange interaction and electron mo- 
bility, thereby shifting the phase boundaries slightly, but 
the overall topology of the ground state phase diagram 
will remain unaffected. Therefore, we believe that at 
least to a first approximation, effects of disorder could 
be neglected. 

Finally, concerning the model (Eq. ^J, there are also 
very interesting issues related to it like the possibility 
of phase separation and canted spin structures, possible 
orbital ordering and low dimensional spin orders which 
have not been investigated here. 



VIII. CONCLUSION 



In conclusion, based on the band structure results 
we have presented an effective model for magnetism in 
Fe doped ZnGa204 which predicts a stable ferromag- 
netic phase when both Fe^"*" and Fe^+ valence states are 
present. If only Fe'^"'" is present -as reported in the Moss- 
bauer spectroscopy^- it is not possible to get ferromag- 
netism via this model, an insulating AFM state would 
have been the most likely ground state. As the system 
is not very dilute and the transition temperature is quite 
high, it is likely that ferromagnetism in this system is 
driven by the kinetic energy of mobile electrons via dou- 
ble exchange rather than interaction between localized 
impurities. A high degree of delocalization of the ex- 
tra electron in Fc^+ could also explain the observation 
of only Fe'^+ states in the Mossbauer experiments. More 
experiments are needed to be done in order to unam- 
biguously detect the states of Fe in ZnGa204 . A study 
of the ground states for a range of doping concentrations 
would be very useful. Photoemission experiments backed 
by detailed first principle calculations are also indispens- 
able to delineate the relevant orbitals that participate in 
the double exchange mechanism. 
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